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We consider a quantum wire with two subbands of spin-polarized electrons in the presence of 
strong interactions. We focus on the quantum phase transition when the second subband starts to 
get filled as a function of gate voltage. Performing a one-loop renormalization group analysis of the 
effective Hamiltonian, we identify the critical fixed-point theory as a conformal field theory having 
an enhanced SU(2) symmetry and central charge 3/2. While the fixed point is Lorentz invariant, the 
effective "speed of light" nevertheless vanishes at low energies due to marginally irrelevant operators 
leading to a diverging critical specific heat coefficient. 

PACS numbers: 71.10.Pm, 64.70.Tg, 75.40.Cx 



As a function of an external gate voltage, the con- 
ductance of quantum wires increases in steps of height 
Go = 2e 2 /h when more and more subbands are subse- 
quently filled. Each conductance plateau transition de- 
fines a quantum critical point where the ground state 
changes its character. Close to each critical point one ex- 
pects universal scaling and a rich behavior as a function 
of gate voltage and temperature. These quantum phase 
transitions and especially their transport signatures are, 
however, only partially understood [H, 0, 0, 0] • 

A quantum phase transition is characterized by the 
critical degrees of freedom, their dynamics, and the un- 
derlying symmetries at the critical point [![. In one- 
dimensional systems, many quantum critical points are 
Lorentz invariant, and conformal field theory (CFT) pro- 
vides a valuable scheme to classify them. In such a case, 
the powerful machinery available for CFTs allows the 
calculation of the critical behavior and the systematic 
investigation of the leading corrections to scaling. In par- 
ticular, a characteristic of CFTs is a linear specific heat, 
C = c(tt/3v)T, inversely proportional to the velocity v of 
excitations, that identifies the central charge c measuring 
the effective number of degrees of freedom [fj, 0] . 

In this Letter we study a model of spinless fermions 
close to the quantum critical point where a second sub- 
band is about to get filled as a function of gate voltage. 
Three different phases @, 0] are relevant to describe this 
system. First, for a small chemical potential, only the 
first subband is occupied, and the electron system forms 
an ordinary Luttinger liquid. Second, for a large chemical 
potential, both subbands are filled, and one obtains two 
Luttinger liquids. There is, however, also a third phase 
which can be visualized as a (fluctuating) zigzag Wigner 
crystal, where the relative motion of charges in the two 
subbands is locked and only one gapless mode survives. 

Upon increasing the chemical potential, the one- 
dimensional phase undergoes a transition to the fluctu- 
ating zigzag Wigner crystal phase 0, 0] • The nature of 



this transition depends on the interaction strength. For 
weak interactions, it has been argued [l], 0] that thermo- 
dynamic properties at the critical point, where the first 
phase becomes unstable, are not qualitatively changed by 
the interactions (see below) . One obtains a Lifshitz tran- 
sition describing the filling of a noninteracting subband of 
fermions with a quadratic dispersion. More precisely, the 
interactions turn out to be dangerously irrelevant. While 
they do not affect the properties at criticality, they nev- 
ertheless induce a zigzag Wigner crystal for small but fi- 
nite filling of the second subband. When the interactions 
are sufficiently strong, however, the universality class of 
the transition changes. The identification of the critical 
properties at such strong interactions is the aim of this 
Letter. 

We show that the corresponding transition has a rather 
unusual structure. On the one hand, we find that the crit- 
ical fixed point is given by a well-known CFT with central 
charge c = 3/2, describing three Majorana fermions with 
equal velocity. As in many other systems the symmetries 
at the critical point are enhanced as neither Lorentz in- 
variance nor the SU(2) symmetry of the c — 3/2 CFT 
are properties of the original model. On the other hand, 
it turns out that, at criticality, the effective velocity — the 
"speed of light" of the CFT — vanishes in the low-energy 
limit. This peculiar behavior arises from marginally ir- 
relevant corrections to the CFT that spoil Lorentz invari- 
ance. It is thus not a property of the fixed point itself, but 
a consequence of the slow flow towards this fixed point. 
We discuss the conditions under which a vanishing of the 
critical velocity occurs generically. A direct consequence 
of this behavior is a divergent critical specific heat coef- 
ficient, 7 = C/T, for T — > 0. 

A known example of a similar situation occurs in a 
model for two-dimensional electrons in graphene interact- 
ing via a long-range Coulomb potential [1, 0] . The Fermi 
points in graphene give rise to an emergent Lorentz in- 
variance that is, however, broken by the static Coulomb 
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interaction. Because of the vanishing density of states 
at the Fermi points, the interaction remains unscreened 
and thus long-range. Although the effective interaction is 
found to be marginally irrelevant, Lorentz invariance re- 
mains broken in the low-energy limit. In contrast to our 
case, the effective velocity is logarithmically enhanced in 
the presence of Coulomb interaction, resulting in a sup- 
pression of the specific heat. 

To be specific, we study the model Hamiltonian 7i = 
H 1 +H 2 +'Hi2, where 



Hi = - I dx 
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The first subband [see Fig. [T{a)] has a finite filling and 
can be treated as a Luttinger liquid [Eq. (TJ]. Its exci- 
tations are plasmons, propagating waves of electron den- 
sity, 8n = ——d x (f), represented by the bosonic field <fi, 
and 9 is its conjugate with [-d x <p(x), 0(x')\ = —iS(x—x'). 
The second subband, on the other hand, contains spin- 
less non-relativistic fermions [Eq. @] that are close to the 
quantum phase transition occurring at r = /i cr — (j, = 
when fj, is increased (e.g., by tuning a gate voltage). Near 
this transition, the interaction among fermions in the sec- 
ond subband is negligible due to the Pauli principle. The 
Galilean invariance of H2 is manifest and implies a crit- 
ical dynamical exponent z = 2 for the fermions in the 
absence of intersubband interactions. However, there are 
two important types of such interactions. The first term 
in Eq. (J3j) is the density-density interaction, g x , and the 
second term, u, arises from pair-tunneling processes of 
electrons as depicted in Fig. Ufa). We assume u > 0. 
Interestingly, the nature of the transition at r = was 
found to depend on the interaction strength Q. 

First, consider the Hamiltonian H without pair tun- 
neling, u = 0. At criticality, r = 0, there is a van Hove 
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FIG. 1: (a) Pair-tunneling process between the two subbands. 
The chemical potential fi is close to the critical value /i cr (r = 
Ha — n)- (b) One- loop RG flow diagram of the effective model 
C e ff determined by the invariant T = X 2 Kvu/(u — v) 4 [see 
Eqs. |[8j| and J9)]. A is the interaction, it is the pair tunneling 
and 11 is the plasmon velocity. There is a line of fixed points 
at it = v parametrized by the dynamical exponent 0. 



singularity in the tunneling density of states (DOS) that 
results in a strong hybridization between the modes of the 
two subbands in the presence of a finite interaction, g x . 
At their band bottom, the fermions are very slow, and 
the fast plasmons can adiabatically adjust themselves to 
the fermions in the second subband. Explicit calcula- 
tions |l[ show that the coupling of the resulting dressed 
fermions to the Luttinger liquid is irrelevant at the criti- 
cal point. Therefore the transition remains a usual one- 
dimensional Lifshitz transition. For fi < /i cr the chemical 
potential is below the band edge, implying a charge gap 
A = /i cr — fi = r. In the limit of a weak interaction g Xl 
this scenario remains valid even in the presence of a finite 
pair tunneling u. However, it was observed in Ref. 0], 
that, for r < 0, the pair tunneling u leads to a ground 
state with a zigzag charge density wave associated with 
a gap, A ~ (—r) a , where the exponent is large, a>l, 
in the weak coupling limit. As a consequence, the actual 
conductance plateau transition does not coincide with 
the Lifshitz transition at r = 0, but is rather shifted to a 
negative value of r < 0, i.e., a finite filling of the second 
subband [|. 

The universality class of the transition at r — changes 
qualitatively when the interaction g x increases. To un- 
derstand the nature of this transition, we apply a uni- 
tary transformation, U — e l f e ^^ dx , to the Hamilto- 
nian such that effectively ip — > U'tpU — ipe~ l6 and 
d x 4> — * dx^+TTip^ip. Switching to a Euclidean Lagrangian 
description, the effective model then can be written as 
£cs = £ll + Asing + Ant, where 
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r is the imaginary time, and A = g x — nv/K. 

For A = 0, the plasmonic and fermionic sectors decou- 
ple, but the fermionic part Asing can now be identified 
with an Ising model that is critical for r = 0. The term 
ip'd x ip in Eq. is negligible close to criticality as it is 
subleading compared with the spatial gradient term in 
Asing- As a consequence, an effective Lorentz invariance 
for the critical fermionic excitations emerges, leading to 
a dynamical exponent z = 1. Note, however, that there 
are two different velocities: the fermions propagate with 
velocity u and the plasmons with velocity v. Moreover, 
the effective model now possesses the self-duality symme- 
try of the Ising model which guarantees that, for A = 0, 
there is a gap A = \r\ in the fermionic sector on cither 
side of the transition. 

In the present work, we analyze the effective model 
Aff perturbatively in the residual interaction A [Ioj |. 
which explicitly breaks the Lorentz invariance of both 
the plasmonic and the fermionic excitations. We show 
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that this Lorentz-invariance breaking leaves traces in the 
critical behavior although the interaction A is found to 
be (marginally) irrelevant in the renormalization group 
(RG) sense. 

Perturbation theory in A is accompanied by logarith- 
mic IR singularities that can be absorbed into a renormal- 
ization group flow of the parameters of the model. Intro- 
ducing a momentum cutoff A, we perform a Wilsonian 
renormalization group step by integrating out pertur- 
batively modes within a momentum shell (A/6, A) with 
b > 1, and we calculate the plasmon and fermion self- 
energy and the vertex correction to one-loop order (see 
Fig. [2|). After rescaling momenta k — > k/b and frequen- 
cies ijj —* w/b z , where z is a so- far-unspecified dynam- 
ical exponent, the momentum cutoff of the perturba- 
tively corrected theory is restored to A. Moreover, in- 
troducing wave function renormalizations, — * y 'Zb 4> 
and ip — * V%F ipi w e can choose the RG conditions such 
that the bosonic velocity v, the Luttinger parameter K, 
and the unit prefactor of ip^d T ip in ((5]) are RG invari- 
ants. This can be achieved by introducing an anomalous 
dynamical exponent 
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A dynamical exponent z / 1 is a manifestation of the 
broken Lorentz invariance. Note that within this RG 
scheme the velocities u and v are not the physical ve- 
locities. In particular, the plasmon velocity w p hy S flows 
under RG (see below). 

The fermionic velocity u, the interaction A, and the ef- 
fective chemical potential r obey the one-loop RG equa- 
tions: 
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These coupled nonlinear differential equations can be 
solved analytically by realizing that the quantity T = 
X 2 Kvu/{u — v) A is an invariant of the RG flow. I deter- 
mines directly the flow diagram shown in Fig. |Tfb). If 
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FIG. 2: One-loop perturbative corrections. The dashed and 
solid lines represent the plasmon and fermion propagator, 
respectively. The dot is the interaction vertex A given in 
Eq. (|6}. (a) Plasmon self-energy, (b) fermion self-energy, and 
(c) vertex-correction. 



the fermionic velocity is larger than the bosonic velocity, 
u > v, the RG flow is towards strong coupling, A — > oo, 
which may indicate either a first order transition or a 
different strong coupling fixed point. If the fermionic 
and bosonic velocities coincide, u — v, the one-loop beta 
functions (JSJ) and © vanish, and one obtains a line of 
fixed points parametrized by the dynamical exponent z 
[Eq. 0]. This line of fixed points is, however, likely 
to be an artifact of the one-loop RG. Fortunately, in a 
two-subband quantum wire we expect the pair-tunneling 
u to be much smaller than the plasmon velocity of the 
first subband v [Icj. In this regime, the RG flow is to- 
wards weak coupling, A — > 0, and is fully controlled by 
our one- loop analysis [ll|. 

If the flow is towards weak coupling, the fermionic ve- 
locity u(b) approaches the bosonic velocity v in the limit 
b — > oo [see Fig. QJb)]. As a consequence, the weak- 
coupling fixed point is characterized by a single velocity. 
Diagonalizing the fermionic sector and refermionizing the 
bosonic mode of the effective model £ e ff for A = and 
U = v, the fixed-point Hamiltonian can be expressed in 
terms of three left- and right-moving Majorana fields, L a 
and R a , respectively: 



^FP = - ^2 V i L ad x L a ~ R a d x R a }- 



(11) 
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At the critical point one has an enhanced SU(2) symme- 
try describin g a rotation in the space of the three real Ma- 
jorana fields 12|. Furthermore, the fixed point is Lorentz 
invariant. Nevertheless, we show below that the residual 
breaking of Lorentz invariance due to the interaction A 
leads to important corrections to the fixed-point Hamil- 
tonian (|lip . although A is formally irrelevant. 

The flow towards the fixed point can be determined 
by considering the limit of a small velocity difference: 
w = 1 — u/w <C 1. By expanding the RG equations ([8]) 
and (El) in lowest order in w. one finds 



dw 
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In this regime, u approaches v very slowly, 
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u(b) w v 1 
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X 3 w. (12) 
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where the prefactor A = (167r 2 /5I) 1 ' 5 has been obtained 
from the exact solution of © and ©. 

If the pair tunneling is very small, u <C v, there is 
an extended regime where the initial flow differs qualita- 
tively from the asymptotic behavior (| 12|1 : 



u(b) 
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An estimate of the crossover scale b* between the two 
regimes yields In b* ~ 1 jX. This scale is very large in 
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FIG. 3: Momentum dependence of the physical velocities of 
the Luttinger liquid, v p i xyB (A/k), and the excitations of the 
second subband, u p h y s(A/A;), for K = 1, |Aj — v = 1 and u = 
0.5 at the quantum critical point. For k — > both velocities 
vanish while u p h ys (A/A;)/i;phys(A/fc) — > 1. Inset: specific heat 
coefficient for various values of u. 



the limit u <C v such that the asymptotic behavior is 
attained only at extremely low energies. 

When interpreting the results given above, it is essen- 
tial to realize that they have been obtained in an RG 
scheme where time and space have been rescaled in an 
asymmetric way as z =/= 1. As velocity is length over 
time, this asymmetric scaling modifies the physical ve- 
locities by an extra factor 



Q(b) = exp 
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with z; p hy S = f2fj itphys = £lu. For low energies, i.e., 
In 6 —>■ oo, z(b) approaches unity because the fixed point 
is Lorentz invariant: 
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However, the flow towards z = 1 is so slow that the 
integral J[z(b') — ljdlnfr' in (| 1 5|) diverges and, thus, 
£1(6 — > oo) = 0. This implies that at the fixed point, 
In 6 = oo, all physical velocities vanish. In Fig. [3]the mo- 
mentum dependence of the two physical velocities it p hy S 
and fphys is shown. As described by Eq. (| 1 3|) the two ve- 
locities approach each other slowly, and both vanish for 
k — ► due to the z renormalization of Eq. (fT5")) . 

The vanishing of the velocities naturally leads to a di- 
vergence of the specific heat coefficient 7. At criticality 
(r = 0), 7 scales as the inverse of the physical velocities, 
"K 6 ) = m(Ej(l + UXb))- Thc RG flow is stopped when 



T(b)/u(b) reaches the cutoff, where temperature obeys 
dT/dlnb = zT. From Eqs. dHJ) and (US), one then ob- 
tains asymptotically 



7 ~ exp 



, 102 , vA 

2 -^5- n — 

2 T 



1/5- 



The inset of Fig. [3] shows the T dependence of 7. 

Away from criticality, r ^ 0, the fermionic excitations 
are gapped. The gap A can be obtained by integrating 
the RG equations (f5|)- (fTU)) up to b = b& with r(b&) — 
m(&a)A. In physical units, this corresponds to the energy 
scale A = m(6a)^(^a)A/6a, where we have again taken 
into account the rescaling of thc time axis. Compared to 
the Ising relation, A = \r\, the resulting T — gap, 
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is suppressed due to thc residual interaction A. The sup- 
pression is, however, slower than power law. The gap A 
should also be observable in the tunneling DOS. 

In this Letter we have shown that the quantum phase 
transition from a Luttinger liquid to a zigzag Wigner 
crystal has a fixed point given by a Lorentz-invariant 
theory of three noninteracting Majorana fermions. How- 
ever, the effective velocities of the modes vanish for low 
energies. This behavior will occur quite generally for sys- 
tems where marginally irrelevant terms break the Lorentz 
invariance. More precisely, whenever the correction to 
the dynamical critical exponent vanishes with 1/ In & or 
slower, the integral in (f!3|) diverges such that — depending 
on the sign of (z — 1) — the effective velocity will either 
diverge or vanish for low energies. 
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